Differential-difference equations associated with the 
fractional Lax operators 

V.E. Adler*, V.V. PostniW 
July 12, 2011 



Abstract. We study integrable hierarchies associated with spectral prob- 
lems of the form Pip = XQip where P, Q are difference operators. The cor- 
responding nonlinear differential-difference equations can be viewed as in- 
homogeneous generalizations of the Bogoyavlensky type lattices. While the 
latter turn into the Korteweg-de Vries equation under the continuous limit, 
the lattices under consideration provide discrete analogs of the Sawada- 

l__l Kotera and Kaup-Kupershmidt equations. The r- matrix formulation and 

^/3 several simplest explicit solutions are presented. 
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(^ 1 Introduction 

r~^| The simplest example studied in this paper is the lattice equation 

«t = u 2 (u 2 ui - u_iu_ 2 ) - u{u\ - u-x) (1) 

where we use the shorthand notations 



(N 



O 



c3 



^ u = u(n,t), u t = dt(u), Uj = u(n + j,t). 

U 

For the first time, this equation was derived by Tsujimoto and Hirota [1, 
eq. (4.12)] as the continuous limit of the reduced discrete BKP hierarchy. 
Recall that both equations 

u,t' = u(ux — U-x) and u t" = u (u2Ux — u-xU-2) (2) 
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are very well known integrable models: respectively, the Volterra lattice 
[2, 3] and the modified Narita-Itoh-Bogoyavlensky lattice of the second or- 
der [4, 5, 6]. One can easily verify that the flows df and df do not commute, 
that is, these equations belong to the different hierarchies. Hence, one should 
not expect a priori that their linear combination remains integrable. Nev- 
ertheless, this is the case: we will show that equation (1) admits the Lax 
representation 

L t = [A,L] 

with the operator L equal to a ratio of two difference operators, namely, 
L = (T 2 + u)~ l {uT 2 + 1)T where T denotes the shift operator Uk —> Ufe+i- 
Equation (1) can be cast into the Hirota's bilinear form which admits a 
family of generalizations depending on a pair of integer parameters (l,m). 
These generalizations were discovered by Hu, Clarkson and Bullough [7, 
eq. (4)] who searched for bilinear equations admitting iV-soliton solutions. 
One of the goals of our paper is to demonstrate that this family of equations 
is associated with the fractional Lax operators of the form 

L = (T m + u)- 1 (uT m + l)T l . (3) 

As usually, any such L is associated with a whole commutative hierarchy 
of equations corresponding to the sequence of difference operators A of in- 
creasing order. We denote this hierarchy dSK^ ,m ', since it can be viewed 
as a discretization of the hierarchy containing the Sawada-Kotera equation 
[8, 9] 

U iT = U 5 + 5UU 3 + 5UiU 2 + 5U 2 U 1 (SK) 

where we denote 

U = U(x,t), U T = d T {U), Uj = di{U). 

For instance, equation (1) belongs to dSK> 1,2 ). The concrete formula of the 
continuous limit in this example is the following, at e — > 0: 



1 £ \rf 

3 + J U ( X - 


4 2e 5 ■ 


- -et,r + ——t 


9 135 . 



u(n,t) = - + —U\x- g et > r + 7^ f J' x = en (4) 

and an analogous formula exists for any (l,m). It should be noted that 
each of equations (2) apart defines a discretization of the Korteweg-de Vries 
(KdV) equation U t t = U^+QUU\ rather than the SK one. Moreover, it is well 
known that actually all Bogoyavlensky type lattices serve as discretizations 
of the KdV equation or its higher symmetries, so that an infinite family of 
discrete hierarchies correspond to just one continuous. Quite analogously, 
the whole family of dSK" ,m ^ hierarchies serve as discrete analogs of the SK 
hierarchy. We hope that this observation makes clear the place of these 
equations in the big picture of integrable systems. 



On the other hand, the differential and difference cases are not quite 
parallel. First, Lax operator for the SK equation 

L = D 3 + UD = (D-f)(D + f)D 

is not fractional. Lax operators given by the ratio of differential operators 
were studied by Krichever [10], however it seems that these examples and 
(3) are unrelated. 

Second, let us consider the problem of discretization for another impor- 
tant example, the Kaup-Kupershmidt equation [11, 12, 13] 

25 

U >T = U 5 + 5UU 3 + —U 1 U 2 + 5U 2 U 1 . (KK) 

Recall that it is associated with the operator 

L = D 3 + UD + hj iX = {D + f)D(D - f) 

and both SK and KK equations are connected through the Miura substitu- 
tions obtained by factorization of Lax operators [14, 15]: 

Usk = f,x - f , Ukk = ~2f,x - f ■ 

Despite of this close relation, it was noted that some properties of the SK 
and KK equations are rather different, see e.g. [16]. It seems that distinc- 
tions between the lattice analogs of these equations are even more deep. A 
discretization of the KK equation is presented in section 4, however, we were 
able to find just one operator L in this case comparing to infinite family (3) 
in the SK case, and no discrete analog of Miura type substitution between 
dSK and dKK is known. 

The contents of the paper is the following. Section 2 contains some 
necessary information on the lattices of Bogoyavlensky type, see also books 
[17, 18]. Section 3 devoted to discretization of the SK equation contains 
the main results of the paper. A general construction of the Lax pairs with 
operator (3) is given in section 3.1. In section 3.3, the r-matrix approach 
in the difference setting [19, 20, 18] is used to obtain explicit formulas for 
the operator A and to prove the commutativity of the dSK' ,m > hierarchy. 
The continuous limit, the bilinear representation, the simplest breather type 
solutions are presented in sections 3.4, 3.5. Section 4 is devoted to discretiza- 
tion of the KK equation and section 5 contains several examples of coupled 
lattice equations associated with more general fractional Lax operators. 



2 Preliminaries 

2.1 Definitions and notations 

We consider differential-difference (lattice) equations of the evolutionary 
form 

u,t = f(um,...,u- m ), u = u(n,t), u )t = d t (u), Uj =u(n + j,t). (5) 

Such equations can be viewed as discrete analogs of continuous evolutionary 
equations like KdV or SK 

U, T = F(U k ,...,U), U = U(x,t), U T = d T (U), U 3 =di(U) 

(the orders m and k may not coincide under the continuous limit). The 
shift operator T : Uj h-> Uj + \ plays the same role for equations (5) as the 
total x-derivative D : Uj i— )■ Uj+\ plays in the continuous case. Differential 
operators are polynomials with respect to D, with the multiplication defined 
by the Leibniz rule DA = D{A) + AD and the conjugation defined by the 
rule D* = —D. In contrast, difference operators are in general Laurent 
polynomials, that is contain powers of both T and T _1 , and the rules for 
the multiplication and the conjugation are TA = T(A)T and T> = T _1 . 
For short, we will use subscripts also for denoting action of T on operators, 
A j = T\A). 

A lattice equation 

V =9(uk,---,u-k) 
is called symmetry of (5) if the compatibility condition Dt(g) = D t i(f) is 
fulfilled, that is 



in 



if,g]* := £ d Us (f)T s (g) - Y, d us(9)T s (f) = 0. (6) 



s=—m 



The lattice is called integrable if it admits an infinite sequence of symmetries 
with the order k greater than any fixed number. The linear space of all 
symmetries is called hierarchy. A conservation law is a relation of the form 

D jt (p(uk, ...,u)) = (T- l)(a(u k +m-i, • • ■ , U- m )) 

which holds true in virtue of equation (5). The discussion of these notions 
and applications to the problem of classification of integrable lattice equa- 
tions can be found in the review article by Yamilov [21]. 

2.2 Bogoyavlensky lattices 

Understanding the structure of dSK(' ,m ' hierarchy is not possible without 
understanding the homogeneous hierarchies of Bogoyavlensky type. A gen- 
eral pattern of (local) equations from dSK"' m ) is given by the formula 

u t = F iL+KM) + ■■■ + F iL+M) + F {L) 



where F^ s > denotes a homogeneous polynomial of degree s with respect to 
the variables Uj and K,L,M are related somehow with the parameters I, m 
and the order k of the flow. Moreover, the first and the last terms in the 
sum always correspond to some (modified) lattices of Bogoyavlensky type 
belonging to the different hierarchies. 

This structure is explained by the following arguments, starting from 
the Lax representation with the operator L (3). Let us consider the scaling 
u — > 5~ m u, T — > ST, then it is easy to see that the limit 5 — > oo sends L 
to the operator L' = u- m T l + T l ~ m and the limit <5 — y leads to L" = 
jrni+i -\- u ~ 1 x. Each of these operators corresponds to its own hierarchy of 
homogeneous lattice equations. The total inhomogeneous equation contains 
both of them together with the intermediate terms which are necessary for 
preserving commutativity of the flows. 

Let us consider the concrete example. One can check that the lattice 

U )t ' = u(wi(Ws + W2 + W\ + w) — W^i(w + W-l + W-2 + W-3) 

- ui(wz + w-i) + u-i(w\ +u>_ 3 )), w := u(l - uiu_i) (7) 

is a higher symmetry of equation (1). Collecting the homogeneous terms 
yields 

U t = F^+ F^ , u t ,= G (7) + G (5) + G (3) 

and the consistency condition of the flows splits to relations 

[F^,G(% = 0, [F^\G(% + [FW,G(% = 0, 

[F (4) i G (3) ]# + [F (2) }G (% = 0, [F( 2 \ G(% = 

where commutator [, ]* is defined by equation (6). As it was already said 
in Introduction, polynomials -F' 4 ^ and F^ 2 > correspond to the modified Bo- 
goyavlensky and Volterra lattices. Polynomials G^ 7 ' and G^ 3 ' correspond to 
their symmetries and the intermediate polynomial G' 5 ^ compensates incon- 
sistency of the hierarchies. 

The Bogoyavlensky hierarchy B' m ) is associated with the operator L = 
T + uT~ m and we recall here several basic formulas regarding this case. A 
detailed theory can be found in the books [6, 18]. More general operators 
of the form L = T l + uT~ m were considered recently in the paper [22]. 

The simplest equation from the B' m ) hierarchy reads 

U t = u(u m -\ h Ul - U-\ U-m). (8) 

This equations and its higher symmetries are associated with the difference 
spectral problem 

and admit the Lax representations 

L ttk = [A^ k \L], L = T + uT- rn , A^ = ir + (L^ m+1 '> k ) (9) 



where n + denotes the projection of any formal series A = ^.- <00 a^'Ti onto 
the linear space of polynomials with respect to T: 

ir + (A)= Y, a ° )ri . 7T-(A) = J2a U) T j . 

0<j<oo j<0 

In particular, 

A (l) = T m+1 + ^ v . = Um + ... + u 

and equation (9) at k = 1 is equivalent to lattice (8). The check is easy: 

L t - [A {1) ,L] = u, t T- m - [T m+1 + v,T + uT- m ) 

= u t T~ m - (u m+1 -u + v- Vl )T - u(v - V- m )T~ m , (10) 

the terms with T cancel and the rest yields the equation. 

In order to prove that equation (9) correctly defines the lattice for any 
k, we have to check that all powers of T except for T~ m vanish in the 
commutator [A(",L]. Since L m+1 is a Laurent polynomial with respect 
to T m+1 , hence A { ' ' is a polynomial with respect to T m+1 . Therefore the 
commutator contains only powers of the form 7 1 ( m + 1 )i+ 1 5 j > _i. On the 
other hand, 

[AM,L] = -[Tr^(L^ m+1 ^),L], 

so that the commutator does not contain positive powers of T and only one 
possible power T~ m remains. 

It can be proven that equations (9) define a special reduction in the Lax 
pair with a generic operator L = T + u^ ' + u^'T^ 1 + • • • + u^ m >T~ m . In 
this case one can choose operators A in the form A = ir + (L k ) with arbitrary 
k. For instance, the Toda lattice hierarchy appears at m = 1. This type of 
multi-field systems was studied, for instance, in papers [20, 22]. 

3 Discretizations of the Sawada— Kotera equation 

3.1 Lax representation 

Let us consider the difference spectral problem 

mp m+ i + 4>i = \(ip m + ml>) (11) 

where m, I are integers. We assume that m, I are positive and coprime, with- 
out loss of generality, since the general case can be obtained by refinement of 
the mesh and/or change of its directions. It is less obvious that the numbers 
m and I can be exchanged: spectral problem (11) is equivalent to 



under the change 

V>(n) = x>(n), A = -x\ [i = -x~ m . (12) 

In the operator form, equation (11) reads 

Pif; = AQ^, P = (uT m + l)T l , Q = T m + u. (13) 

The isospectral deformations are defined by equation ipj = Aip with some 
difference operator A. The corresponding Lax equation 

L t = [A,L], L = Q- 1 P (14) 

can be rewritten as the system 

P tt = BP- PA, Q )t = BQ-QA (15) 

where one of equations can be considered just as a definition of B. Let P, Q 
be as in (13), then this system is equivalent to equations 

u t = B(T m + u) - (T m + u)A, 
B(T 2m - 1) = A m T 2m - Ai + uAT m - uA mM T m . (16) 

In order to resolve the latter we make the assumption that operator A is of 
the form 

A = F{T m - T~ m ) (17) 

then B is found as the difference operator 

B = F m T m -F l T- m + u{F-F m+l ) (18) 

while first equation (16) turns into 

u t = T m Fu + uFT- m -uT m F t - F l T~ m u + F m - F t + u(F - F m+l )u. (19) 

It is clear that the same evolution of the variable u is defined by the conju- 
gated operator F^ and, moreover, all terms T J ', j \ m can be thrown away. 
This means that we can find F as a self-adjoint operator F = F> which is a 
Laurent polynomial with respect to the powers T m : 

F = f( k )T km -\ |_ f( 1 )j im _|_ f(°)_)_7 1 - m f( 1 )_| ^_rp~kmr{k) £ > Q ^20) 

Certainly, the coefficients depend on k,l,m, so that it would be more rig- 
orous to write fv' k ' l > m ) instead of /"•', but we will consider these numbers 
fixed at the moment. 

Collecting the coefficients at T jm , j > 0, yields the relations 

j-i) _ „ ftf-i) = f.(?) _ f(i) + ,,.,.,_ (f& . - fUh 



1 jmJ m u J m +l ~ H Jm ^ au jm{J m+ i J 

I- f { 



+ u jm 4 j+1) -uf^ +1 \ j = l,...,k + l, (21) 



where it is assumed for convenience that /«) = at j > k. The coefficient 
at T° gives an evolutionary equation for u: 

u t = 2u(/« - /«) + n 2 (/(°) - fMj + /<?> - /<°>. (22) 

System of equations (21), (22) defines the k-th. flow in the hierarchy dSK> l ' m >. 
If we are interested in the local evolution only then we require that all 
fU' can be recurrently found as functions of an finite set of variables Ui. In 
this case a certain restriction on the values of k appears and a part of the 
flows is rejected. Indeed, consider equation (21) at j = k + 1, 

u {k+1)m ft ] = uf%, (23) 

or 

(T l - l)(log/W) = {T^ m - l)(logu). 

It can be proven that it is solvable with respect to f^ ' if and only if (k-\- \)m 
is divisible by I and the solution is, up to a constant factor, 

f {k) = u- m ui_ m ---u {s _ 1)l _ m , (k + l)m = si. (24) 

Since / and m are coprimes, hence the local flows may appear only if k = 
pi — 1 and s = mp. The fact that the rest equations (21) for such k are 
solvable indeed will be verified later in section 3.3. The case I = 1 is the 
only one when there are no restrictions on k and the simplest choice k = 
brings in this case to the following family of lattices. 

Theorem 1. For any m > 0, the simplest equation in the hierarchy dSK\ l,m ' 

u t t = u 2 (u m ■••ui — u_i • • • u_ m ) - u(n m _i •••«! — U-l • • • U\-m) (25) 

possesses Lax representation (14) with the operators 

p = u T m+1 + T, Q = T m + u, 
A = f(T~ m -T m ), B = hT- m - f m T m + u(f m+1 - f) 

where f = U—\ ■ ■ ■ u_ m . 

Proof. A direct computation (cf with (10)) proves that both equations (15) 
with given P,Q,A,B are equivalent to relations 

Umfm = Ufm+lj u ,t = U (fm+1 ~ f) ~ fm + /l- 

The former defines the variable / (up to a constant factor) and the latter is 
equivalent to lattice (25). □ 



In particular, equation (25) at m = 2 coincide with (1) and at m = 1 it 
is just the modified Volterra lattice 

u )t = u 2 (m - u-i). 

It should be remarked that gauge equivalence (12) between the spectral 
problems can be extended on the level of nonlinear equations and the same 
flow (25) appears also as a member of dSK( m ' 1 ' hierarchy. However, op- 
erator (20) is much more complicated in this case: it contains all powers 
T" 1 - 1 , T m ~ 2 , ..., T l ~ m comparing with just F = /(°) in dSK^ 1 '" 1 ) case. 

Computing of higher symmetries quickly becomes involved, because find- 
ing of F requires (discrete) integration of rather bulky expressions. For in- 
stance, the second flow in the hierarchy dSK' 1 '" 1 ) is, according to (22), of 
the form 

V = 2u(/« _ /«) + tt * (/ (0) _ /W i} + f io) _ A (0) 

where functions f^\f(°' are defined by relations 

Uo f (i)_ uf (i) u f (0)_ f (o) _ -(l) _ -(l) _ ff (i) _ f (i) s 

a 2mJ m — u Jm+l> u rnJ m u Jm+l ~ Jl Jm uu m\J Jm+lJ- 

This yields, up to integration constants, 

/ (1) = M m _i • • • li_ m , / (0) = (w-\ h W- 2m +l)U-l ■ ■ ■ U^ m , 

w := (1 - u m -iU-.i)u m -.2 ■■■uo 

(at m = 2 equation (7) appears). One can check straightforwardly that the 
obtained flow commutes with (25) indeed. A general proof and a way to 
bypass the integration are given below in section 3.3. 

Adopting nonlocal variables leads to some extension of the hierarchy. In 
this case we consider equation (23) as a constraint which defines the variable 
j{ k ) f or an y fa. Then we arrive to the following system which generalizes (25) 
for any i, making the picture more uniform. We will return to this system 
in section 3.5. 

Theorem 2. For any coprime m,l, the simplest system in the extended 
dSKS l ' m ^ hierarchy 

u m f m = uf m+ i, u yt = u 2 (f-f m+ i) + f m -fi (26) 

possesses Lax representation (14) with operators 

P = uT m+l + T l , Q = T m + u, 
A = f(T- m - T m ), B = hT- m - f m T m + u(f m+l - /). 



m = 2 : u t t = u {u2U\ — U-1U-2) — u{u\ — U—i) 

U = V\V V } t = V\V^V-\(v2Vl — f_l1>_2) — V 2 (V\ — V-l) 

m = 3 : u t t = u (U3U2U1 — u-iU-2 u -3) — u(u2Ui — u_iu_2) 

V = U\U Vj = v(vsVl + V2V — VV-2 — V-lV-z) — V(V2 + V\ — V-\ — V-2) 

U = V2V\V Vj = V2V 1 V V_ 1 V-2(vsV2Vl — V-IV-2V-3) 

— -uif 3 -u_i(i;2fi — V-xV-2) 

m = 4 : U( = U 2 (u4UsU2U\ — M_iU_2'"-3^-4) — u(u-^U2U\ — -U_iU_2^-3) 

U = V2V V t t = V2V\V V-\V-2{V4V$V2V\ — V -\V -iV --&V -4) 

- V\V 2 V-\{v?,V2V\ - V-1V-2V-3) 

Table 1. Examples of lattices (25) from dSK^ 1 '" 1 ) and their modifications 

3.2 Modified lattices 

Equations under consideration can be rewritten in several ways by use of 
difference substitutions. The simplest kind of substitution is introducing a 
potential. Let A be a constant operator, then substitution u = A{v) maps 
solutions of equation v t t = f[A(v)] into solutions of equation u t t = A(f[u]). 
Table 1 contains several instances of such kind, up to the change u — > e u , 
v — > e v . 

Another kind of substitutions are Miura type transformations. Let <p be 
a particular solution of spectral problem (11) corresponding to a value X = a 
of the spectral parameter. Then one readily finds that the ratio h = pi/<p 
is related with the potential u by formula 

ah m -i---h-hi-i---h 

M : u= . 

hm+l-i ■■■ n- a 

This defines a difference substitution, according to the following statement. 

Theorem 3. Let u satisfies an equation (22) from dSKS ' m > , then h = ip±/ip 
also satisfies a lattice equation which can be written as a conservation law 

(logh),t = (T-l)S[h\. (27) 

Proof. Since p is governed by equation <p t t = Ap = F((p m — p- m ), hence 

(logh), t = {T- l)(log^), t = (T- l)(-F{ip m - v?_„ 
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Coefficients of the operator F are functions on the variables hj , being func- 
tions on Uj's. The ratios of the form tpk/<P can be expressed through hj as 
well and therefore an equation of the form (27) holds. □ 

It is worth noticing that an infinite sequence of conservation laws for the 
original lattice (22) can be obtained from (27) by use of the classical trick 
with the inversion of Miura map u = M~(h, a) as a formal power series with 
respect to a [23]. 

Second Miura map is obtained by replacing h — > 1/h, a —> 1/a which 
results in the mapping 

, r + a/i m _M_i • • • hi — hm+i-i ■■■h m 
M T : u = 



hm+i-i -"h-a 

This substitution relates the same equations as M~, due to invariance of 
the spectral problem with respect to the change n — > —n, X — > 1/A. There- 
fore, the composition M~(M + )~ 1 defines a Backhand transformation which 
relates two copies of the dSK"' m ) hierarchy. Recall that Backhand transfor- 
mation for the continuous SK equation was derived in [24]. 

A particular example at I = 2, m = 1 is given by substitutions 

M-= u= [°'-' ll)h , M+: u= h J 2 [ a - hl) 
h2h\h — a hih\h — a 

which map solutions of the modified equation 

h(a — h) ( h(a — h{)(a, — h-\)(h2h\ — h-\h-2) 
h\hh—\ — a\ {hih\h — a)(hh—ih—2 — a) 

into solutions of (1). 

3.3 r-matrix formulation 

In this section we prove that: 

(i) if the constraint (23) is resolved by formula (24) then the further 
recurrent relations (21) are solved in the local form as well, so that the 
(local) hierarchy dSK'' ,m ^ is correctly defined; 

(ii) the flows corresponding to the different k commute. 

In achieving this goal the r-matrix approach is an indispensable tool, see 
e.g. [19, 20, 18]. Let us consider the Lie algebra of the formal Laurent series 
with respect to the powers T m of the shift operator: 

g (m) = (^ g (J) T Jm"\ 

j<oo 

with the commutator [A, B] = AB — BA. It is easy to see that any element 

q = (fc+1) ji(fe+l)m _|_ (k)j,km _|_ (fc-1) jr(fc-l)m _)___ 
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of this Lie algebra admits an unique decomposition of the form 

G = F(T m - T- m ) + H (28) 

where F = F' is a self-conjugated difference operator and H is a formal 
series which contains only nonpositive powers of T m . Each of the linear 
spaces 

g M = | F(Tm _ r - m) | p = F t} ; fl (m) = | J- ^(jj^m j 

i<o 

constitutes a Lie algebra: for g_ this is obvious and for q™ we have 

\F(T m — T~ m ) p'(p m — T~ m )] = ( P + P^](T m — p~ m ) 

where P = F(T m - T- m )F'. 

Thus, formula (28) is the decomposition (in the vector space sense) 

M = gM e L m) 

of the Lie algebra into the direct sum of two Lie subalgebras. This decom- 
position defines the projections n± on the g^™ component and the r-matrix 
r = 5( 7r + — 7T— )■ Now we can formulate the following theorem about Lax 
equations (13), (14) with fractional L operator. 

Theorem 4. Let l,m be coprime, P = (uT m + l)T l , Q = T m + u and let 
L = Q~ 1 P be expanded as a formal Laurent series. Then the flows 

L tp = [n + (If m ),L] (29) 

are correctly defined for all p = 1,2, ... , coincide with the dSK^ ' m > flows 
introduced by equations (21), (22) and commute with each other. 

Proof. After expanding, L takes the form 

L = (1 - n_ m T" m + (u. m T- m f -... )(u- m + T~ m )T l 



U- 



x T l + (1 - u. m u. 2m )T l - m + ... . 



Differentiating this series turns (29) into an infinite system of equations for 
a single variable u, and the correctness means that all these equations must 
coincide. To prove this, we compare representation (29) with Lax equation 
(14) in fractional form. 

Notice that L itself does not belong to the Lie algebra q^-" 1 ' , but its power 
G = L pm does, so that the projection A = vr + (G) = F(T m - T~ m ) makes 
sense. We denote the order of operator F as k = pi — 1, in agreement with 
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(20) and (24). The coefficients of F are uniquely computed from coefficients 
of G accordingly to the recurrent relations 

/ (fc+2) = / (fc+1) = 0, /0) =5 (i+i) + / (i+2) ) j = k ,k-l,...,0 

so that all coefficients are local functions of Uj (in particular, f( k > is given 
by (24)). Moreover, the order of (29) right hand side is equal to I, because 
[n+(G),L] = — [tt-(G),L\. This proves that F provides a solution of the 
recurrent relations (21) as well (which is unique up to integration constants). 
Indeed, these relations were derived from the condition that terms with 
T' +1 ' m ,. . . ,T m in equation (19) cancel which is equivalent to cancellation 
of the powers j 1 ( fc + 1 ) m +' )- . . ^pn+l - m fo e original Lax equation (14). Thus, 
flow (29) coincides with a flow from dSK' ' m ' which is, therefore, local. On 
the other hand, this proves correctness of (29), since the whole infinite set 
of equations turns out to be equivalent to the single equation (22). 

The proof of the commutativity is standard. Let G = L pm and A! = 
7T+(G") then 

(L tp ),t p , ~ (L, tp/ ),t p = [\, - A' tp + [A,A'],L], 

so it is sufficient to prove that 

Atpl -A' tp + [A,A'] = 0. 

Since At , = 7r+(L4', G]) and [G, G'\ = 0, this is equivalent to 

7C+([A',G]-[A,G'] + [A,A'}) 

= %+([G' - tt-{G'), G]-[G- tt_(G), G'\ + [G- tt_(G) ; G' - %-(G')}) 
= 7r+([7r-(G),7r_(G')])=0 

as required. □ 

3.4 Continuous limit 

Here we compute the continuous limit for the basic flow of the extended 
hierarchy dSK"' m ) defined by equation (26). There is a certain technical 
difficulty in the prolongation of the continuous limit on the variable / which 
is not local at I ^ 1. In order to solve the constraint, this variable should 
be considered as a series with respect to the small parameter. Up to this 
complication the continuous limit is very similar to example (4) from Intro- 
duction. We postulate that, at e — > 0, the variables u, f are of the form 

u(n, t) = a + abe 2 U(x + cet, r + de 2 t), 

oo 

f(n,t) = l + ^2e s Y s (x + cet,T + de 2 t), x = en (30) 



s=2 
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with undetermined coefficients a, b, c, d. Functions Y s are expressed through 
the function U and its partial derivatives with respect to x after substituting 
into first equation (26) and taking the Taylor expansion about e = (clearly, 
one can neglect the dependence on t here). We find, omitting the unessential 
integration constants: 

Y 2 = ft/, 

m(m + l)b 
13 = y U lt 

m(m + l)(2m + l)b rr m(m — l)b 2 o 

n = -* — a l U2 + -^^t/ 2 , 

m 2 (m + l) 2 b m(m 2 — l 2 )b 2 

Yb = 24/ U3 2P UUl > 

m(m + i)(2m + i)(3m 2 +3ml-l 2 )b m(m 2 - l 2 )(3m + 2l)b 2 2 

Y6 = 7201 U * + 24f 2 ^ 

m(m 2 — l 2 )(2m + l)b 2 TTTT m(m — l)(m — 2l)b 3 r ■> 

+ ^ ^ } -uu 2 + -J ^ L f/ 3 

This is enough, since we need only terms up to e 7 when substituting into 
second equation (26). The coefficients a,c are found from the requirement 
that the low order terms vanish while the coefficients b, d are responsible for 
the scaling of U and t and can be chosen arbitrarily. Finally, we come to 
the following statement. 

Theorem 5. Continuous limit (30) with the values of parameters 

m — I , ml , m 3 (/ 2 — m 2 ) 

— • — c = 2m, a - 



m + V 6 180 

sends systems (26) into the Sawada-Kotera equation 

U T = U 5 + 5UU 3 + W 1 U 2 + 5U 2 Ul 

The higher flows of the SK hierarchy can be derived analogously from 
suitable linear combinations of the dSK^' m ) flows. However, the general 
formulas become rather complicated and we restrict ourselves by the fol- 
lowing concrete example corresponding to the local hierarchy dSK' 1 ' 2 ). Let 
u t t = 88u t t 1 + 27u t t 2 where the flows dt t and dt 2 are defined by equations (1) 
and (7) respectively, then the formula 

/ n 1 £2 tt( 200 16e 7 \ 

u(n, t = — | U x et, t 1 , x = en 

v ' ; 3 9 V 9 ' 189 / 

defines the continuous limit to the 7-th order symmetry of SK equation 

U T = U 7 + 7UU 5 + 14t7iZ7 4 + 21U 2 U 3 + UU 2 U 3 

28 
+ 42UU 1 U 2 + 7Uf + —U 3 U 1 . 

o 
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It is well known that there are gaps in the sequence of orders k of equations 
from the SK hierarchy, namely, the restrictions k \ 2, 3 are fulfilled, so that 
the next higher symmetry is of 11-th order. The natural question appears, 
how this agrees with relations (20)-(22) or (29) which show that in the 
discrete case there are no gaps multiple 3. It turns out that their appearance 
is an artefact of the continuous limit. A straightforward computation shows 
that if we consider a linear combination with the next dSKA 1,2 ' flow u t t = 
u^j + au 7 t 2 + Pu t t 3 and set 

u(n, t) = a + be 2 U(x + cet, r + de 9 t), x = en 

then all parameters are uniquely determined by the condition of vanishing 
the terms up to e 10 , however then the coefficients at e 11 cancel automatically 
and only the trivial flow U T = appears. 

3.5 Bilinear equations 

The constraint (23) can be solved by introducing additional variables and 
this leads to a convenient representation of the basic system (26) of the 
extended dSK^ ,m ) hierarchy. Let 

V[ v 

u= —, f = 

V V- m 

then first equation (26) is satisfied identically and the second one is equiva- 
lent to 

(T*-1)^ = (T™-1)(^ V -L-). 

V \Vl- m V- m ) 

Further substitutions 

w m w m+ iw w m w- m 

v = =>- u= , / = k — 

w w m wi w z 

bring to the bilinear equation 

w l} tw-wiw,t = w m wi- m -w- m wi +m . (31) 

For the first time, it appeared in paper [7], in a slightly more general form 

wijw - wiw t t = w m wi- m - aw- m wi +m + fiwwi 

which is reduced to (31) by the point change w(n,i) = e^*a n w(n,t). In 
particular, it was proven in [7] that this equation admits iV-soliton solutions. 
Here, we consider in more details a specification of 2-soliton formula which 
leads to the breather solution. 

The substitution of the 2-soliton Ansatz 

w(n,t) = 1 + ei + e 2 + Ai 2 eie 2 , e« = g™exp(-Wjt + <5j) (32) 
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Figure 1. The values of q — pe lip inside the bounded domains in C 

correspond to the regular potentials u(n,t). The values along the dashed 

lines correspond to the potentials periodic in t. 



into (31) gives us the dispersion relation and the phase shift: 

(J -q l )(q m -q m ) 

The direct check proves that then the 3-soliton Ansatz 

w = 1 + e\ + €2 + e% + A\ 2 e\e2 + Ai 3 eie3 + A 23 e2e3 + Ai 2 A 13 A23eie2es 

satisfies (31) automatically. It is interesting to compare these formulas with 
their counterparts for the continuous SK equation [8, 9, 25, 26] 

\2/„2 



(Xj - Xj) J (x/ - XjX,' + x^ 



ei = exp(xjx - Wjt + Si), uji = xf , A i:j - ,.,, .-, — 

(xj + >CjY\ycf + XiXj + x~) 
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Figure 2. A moving and a stable breathers. The values of parameters: 
p = 1.2, if = 2ir/3 (left); p = 1.6, p = 37r/4 (right); in both cases I = 1, 

m = 2, a = P = 0. 



Formula (32) allows us to obtain the breather-type solutions as well, if we 
choose 



9i =pe 



V 



92 = pe 



-iv 



a + i/3, 



a — i/3. 



The regularity of the potential u(n, t) is achieved under certain restrictions 
on the value of q. In order to show this, rewrite relations (33) as follows: 

oo = p + w, [i = [p —p )cosrrnp, v = (p +p jsmrrup, 

Ap m+ sin lip sin imp 



A 
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(1 - p 2l ){l - p 2m ) ' 

then a simple algebra brings (32) to the form 



w 



1 + 2z cos y + A12Z , y = ipn — vt + /3, 



p n e a -^. 



In particular, if <p = 2 ^ 7r then p = and solution u is periodic in t. The 
necessary and sufficient condition for u to be regular is that the function 
w does not vanish at any n, t. In the generic case the variables y, z are 
independent and then this is equivalent to the condition that the trinomial 
1 + 2z + A\2Z 2 does not vanish at real z, that is 

(p l - p~ l )(p m ~ p' m ) +4 sin lip sin imp < 0. 

Thus, we see that already two-phase solutions in these models exhibit a 
nontrivial zone structure of the spectrum. The corresponding domains in 
the plane q = pe lip are shown on fig. 1, and the examples of solutions u(n, t) 
are shown on fig. 2. 
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4 A discrete analog of the Kaup— Kupershmidt equa- 
tion 

The Kaup-Kupershmidt equation 

U T = U 5 + 5UU 3 + y£/i£/ 2 + 5t/ 2 C/i 

is associated with the spectral problem Lip = Xif) where L is the skew- 
symmetric ordinary differential operator of third order 

L = D 3 + UD + hj tX = {D- f)D(D + /), U = 2f jX - f. 

When we find a discrete analog, a difficulty is that a symmetric or skew- 
symmetric difference operator can be of even order only. A way to overcome 
this is to consider a 6th order difference problem, but on the odd nodes of the 
lattice only, so that effectively it is of 3rd order with respect to the double 
shift T 2 (however, the coefficients may depend on the variables associated 
with the even nodes as well). Let us consider the spectral problem 

u_ 3 V>-3 + V'-i = A(^i + uip 3 ) (34) 

or, in the operator form, denoting K = uT 3 + T: 

iTty = XKip. 

The Lax equation for the operator L = K~ l K^ can be written in the form 
of system (15). It admits the reduction B = —A' which yields the equation 

K tt + A ] K + KA = 0. (35) 

The operator A is found as a Laurent polynomial with respect to the even 
powers of T, 

A = a {k) T 2k + ■■■ + a(-^T- 2k , 

and a direct analysis of equation (35) at k = 1,2 proves the following state- 
ment. 

Theorem 6. Equation (35) with K = uT 3 + T is equivalent to the nonlocal 
lattice equation 

u,t! = u{f 2 u 2 - /i«i + f-iu-i - f-2U-2) + h ~ f-i, hu = f-\u 2 

under the choice 

A = -fT 2 + jl 2 u_ 2 - /_!«_! + f- 3 T- 2 ; 

and it is equivalent to the local lattice equation 

U t t 2 = u(v 3 - V2 + Vi - V-i + V-2 - V-3 - U 2 + U- 2 ) , V := U\UU-i (36) 
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under the choice 

A = uiT A - u-iT- 4 + (1 - u-iu~ 2 ){T 2 - T~ 2 ) 

+U-\ — U-2 — V + V-l — V-2 + V-3. 

It is worth noticing that, alternatively, one can use the following pair of 
operators (cf with the gauge equivalence (12)): 

K = uT 3 + T~ 1 , A = -uiii_iT 4 + u_ 2 n_ 4 T- 4 -t; + v_i-U-2 + V-3- (37) 

The continuous limit to the KK equation is of the same general form as 
before, namely, for the flow (36) it reads 

x 1 4 2tt( 8 64e5 \ 

u{n, t 2 ) = - + -e U[x- -et 2 , r + j^rte J , x = en. 

5 Examples related to generic operators 

Recall that, according to [20], the Bogoyavlensky type lattices can be viewed 
as reductions of more general multi-field models associated with the spec- 
tral problems Lijj = A^ for generic difference operators L = u^ m 'T m + 
u (m-i) T m-i H 1_ u (i-i) T i-i + u (-i) T ~i_ Here m j are any p OS itive in- 
tegers, and one can adopt the normalization w m > = 1 or u"' = 1 without 
loss of generality. A part of the flows from the corresponding hierarchy is 
consistent with the constraints u^ m ~ 1 ' = • • • = vS l ~ l > = and this reduc- 
tion brings to the Bogoyavlensky lattices. A detailed study of some other 
reductions can be found in [22]. 

The lattices introduced in the previous sections are related with the 
spectral problems Pip = XQvp where operators P, Q are binomial. It is 
natural to expect that these lattices also define reductions for some multi- 
field equations related with more general operators P, Q. The study of such 
models is beyond the scope of the present paper and we restrict ourselves 
by three typical examples. 

Example 1. First, let us consider the Lax equations Pj = BP — PA, Qj = 
BQ — QA for the binomial operators P, Q with different potentials: 

P = uT 3 + T, Q = T 2 + v. 

If v = u then operators A, B are given by formulas (18), (20) with a self- 
adjoint operator F which contains only even powers of T 2 . In the general 
case two sets of operators A, B appear, containing positive or negative pow- 
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ers of T 2 . The simplest operators and corresponding flows are the following: 

A- = v- 2 v-!T- 2 + /_ 3 + f- 2 , B~ = v^vT- 2 + /_l + /, 

u,t- = u{f-i - fl), 

v ,t- = V(f + /-l - f-2 - f-3 ~Vl + V-l), f := uviv 2 ; 

A + = U-2U-1T 2 + £_i + g, B + = uuiT 2 +g + gi, 
u }t + = u(g + 51 - 52 - 93 - u-l + wi), 
v,t+ =v(9i -9-1), g:=u- 2 u-iv. 

The flows d t - and d t + commute, and the flow dt = d t - — d t + admits the 
reduction v = u which brings to the dSK equation (f). It should be noted 
that the same flows can be obtained starting from the gauge equivalent 
operators P = uT 3 +T 2 , Q = T + v. 

Example 2. Now let us consider trinomial operators 

P = U T 3 + pT 2 + T, Q = T 2 + qT + v. 

In this case operators A, B contain the odd powers of T as well. The simplest 
operators and the corresponding flows are of the form 

A~ = V-\T~ X + V-1P-2, B~ = vT^ 1 + v\p, 
u >x - = u(u-iq - uiq 2 -p + Pi), 
P,x- = P(u-1Q ~ uqi) +u-U-!, 
v >x - = v(u-iq - u- 2 q-i), 
q x - = uv\ — u- 2 v; 

A + = u- 2 T + u- 2 q-i, B + = uT + U-\q, 
u >x + = u(v 1 p-v 2 pi), 

P jX + = U-\V - uv 2 , 

v >x + = v(vip - V-ip- 2 + q-i - q), 

Q,x+ = q(viP ~ vp-i) +v-vi. 

Example 3. Let us consider the following generalization of the spectral prob- 
lem (37): 

K^ip = XKip, K = uT 3 + v- 1 T + T~ 1 . 

The isospectral deformations are defined by the operators A = a^'T + 
a (k-l)rp2k-2 _|_ . . . _|_ a(~ k ">T~ 2k . The simplest case k = 1 results in 

A = U-\T — ii- 2 T~ + U-1V-1 — U- 2 V- 2 
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and equation K^ + A^K + KA = is equivalent to the lattice 

U yt = -u(u 2 V2 - UiVi + U-lV-l - U- 2 V-2 — Vl + V-l), 
V t t = —v(u\Vi — U-lV-l) + U2U1 — U-1U-2 + U\ — U-\. 

The higher symmetry corresponding to k = 2 is too bulky and we do not 
write it down, however one can check that it admits the reduction v = to 
the dKK equation (36). In contrast, the flow dt itself does not admit this 
reduction. 

6 Conclusion 

In this article we introduced a family of integrable lattice hierarchies associ- 
ated with fractional Lax operators. In particular, these hierarchies contain 
equations found earlier in [1, 7] by use of the Hirota bilinear formalism. 
We proved that these equations serve as semi-discrete analogs of SK and 
KK equations. An important question which remains open is about the 
Hamiltonian structure of the presented equations. As usually, the existence 
of Lax representation allows to obtain a set of conserved quantities which 
presumably are Hamiltonians, and moreover, the applicability of r-matrix 
approach suggests that some more or less standard Poisson bracket should 
exist. However, no explicit answer is found yet. Another intriguing question 
is about possible relations with the models introduced in [27, 28] within the 
theory of the lattice W algebras. 

Acknowledgements 

We are grateful to Ya.P. Pugay, Yu.B. Suris and A.K. Svinin for many 
stimulating discussions. The research of V.A. was supported by grant NSh- 
6501.2010.2. 

References 

[1] S. Tsujimoto, R. Hirota. Pfaffian representation of solutions to the dis- 
crete BKP hierarchy in bilinear form. J. Phys. Soc. Japan 65 (1996) 
2797-2806. 

[2] V.E. Zakharov, S.L. Musher, A.M. Rubenchik. Nonlinear stage of para- 
metric wave excitation in a plasma. JETP Lett. 19:5 (1974) 151-152. 

[3] S.V. Manakov. Complete integrability and stochastization of discrete 
dynamical systems. Soviet Physics JETP 40:2 (1975) 269-274. 

[4] K. Narita. Soliton solution to extended Volterra equation. J. Phys. Soc. 
Japan 51:5 (1982) 1682-1685. 

21 



[5] Y. Itoh. Integrals of a Lotka-Volterra system of odd number of variables. 
Progr. Theor. Phys. 18 (1987) 507-510. 

[6] O.I. Bogoyavlensky. Algebraic constructions of integrable dynamical 
systems — extensions of the Volterra system. Russian Math. Surveys 
46:3 (1991) 1-64. 

[7] X.B. Hu, P.A. Clarkson, R. Bullough. New integrable differential- 
difference systems. J. Phys. A 30:20 (1997) L669-676. 

[8] K. Sawada, T. Kotera. A method for finding n-soliton solutions of the 
KdV equation and KdV-like equations. Progr. Theor. Phys. 51:5 (1974) 
1355-1367. 

[9] P.J. Caudrey, R.K. Dodd, J.D. Gibbon. A new hierarchy of Korteweg- 
de Vries equations. Proc. Roy. Soc. London A 351 (1976) 407-422. 

[10] I. Krichever. Linear operators with self-consistent coefficients and ra- 
tional reductions of KP hierarchy. Phys. D 87:1 4 (1995) 14-19. 

[11] D.J. Kaup. On the inverse scattering problem for cubic eigenvalue prob- 
lems of the class tp xxx + 6Q ip x + dRip = \ip. Stud. Appl. Math. 62 (1980) 
189-216. 

[12] B.A. Kupershmidt. A super KdV equation: an integrable system. Phys. 
Lett. A 102:5 6 (1984) 213-215. 

[13] V. G. Drinfel'd, V. V. Sokolov. Lie algebras and equations of Korteweg- 
de Vries type. J. of Math. Sciences 30:2 (1985) 1975-2036. 

[14] V.V. Sokolov, A.B. Shabat. (L, j4)-pairs and a Ricatti type substitution. 
Fund. Anal. Appl. 14:2 (1980) 148-150. 

[15] A. P. Fordy, J.D. Gibbons. Some remarkable nonlinear transformations. 
Phys. Lett. A 75:5 (1980) 325. 

[16] M. Musette, C. Verhoeven. Nonlinear superposition formula for the 
Kaup-Kupershmidt partial differential equation. Physica D 144 : 1>£ 
(2000) 211-220. 

[17] O.I. Bogoyavlensky. Breaking solitons. Nonlinear integrable equations. 
Moscow: Nauka, 1991. 

[18] Yu.B. Suris. The problem of integrable discretization: Hamiltonian ap- 
proach. Basel: Birkhauser, 2003. 

[19] A.G. Reyman, M.A. Semenov-Tian-Shansky. Integrable systems. 
Group-theoretical approach. Izhevsk, 2003. (in Russian) 



22 



[20] M. Blaszak, K. Marciniak. i?-matrix approach to lattice integrable sys- 
tems. J. Math. Phys. 35:9 (1994) 4661-4682. 

[21] R.I. Yamilov. Symmetries as integrability criteria for differential differ- 
ence equations. J. Phys. A 39:45 (2006) R541-623. 

[22] A.K. Svinin. On some class of homogeneous polynomials and explicit 
form of integrable hierarchies of differential-difference equations. J. 
Phys. A 44 (2011) 165206 (16pp). 

[23] R.M. Miura, C.S. Gardner, M.D. Kruskal. Korteveg-de Vries equation 
and generalizations. II. Existence of conservation laws and constants of 
motion. J. Math. Phys. 9:8 (1968) 1204-1209. 

[24] J. Satsuma, D.J. Kaup. A Backlund transformation for a higher order 
Korteweg-de Vries equation. J. Phys. Soc. Japan 43 (1977) 692-697. 

[25] E. Date, M. Jimbo, M. Kashiwara, T. Miwa. KP hierarchies of orthog- 
onal and symplectic type. Transformation groups for soliton equations. 
VI. J. Phys. Soc. Japan 50 (1981) 3813-3818. 

[26] A. Parker. A reformulation of the 'dressing method' for the Sawada- 
Kotera equation. Inverse Problems 17:4 (2001) 885-895. 

[27] Ya.P. Pugay. Lattice W algebras and quantum groups. Theor. Math. 
Phys. 100 (1994) 900-911. 

[28] K. Hikami. Generalized lattice KdV type equation reduction of the 
lattice W 3 algebra. J. Phys. Soc. Japan 68 (1999) 46-50. 



23 



